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Abstract 

I provide evidence that the 2D RP N ^ model for N > 3 is 
equivalent to the 0(iV)-invariant non- linear cr-model in the con- 
tinuum limit. To this end, I mainly study particular versions of 
the models, to be called constraint models. I prove that the con- 
straint RP N ~ 1 and O(N) models are equivalent for sufficiently 
weak coupling. Numerical results for the step-scaling function of 
the running coupling g 2 = m(L)L are presented. The data con- 
firm that the constraint 0(N) model is in the same universality 
class as the O(N) model with standard action. I show that 
the differences in the finite size scaling curves of RP 1 ^^ 1 and 
O(N) models observed by Caracciolo et al. can be explained as 
a boundary effect. It is concluded, in contrast to Caracciolo et 
al., that RP N ~ l and O(N) models share a unique universality 
class. 
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1 Introduction 



Motivated by the close analogy with non-Abelian lattice gauge theories in 
four dimensions non-linear cr-models in two dimensions have been studied 
intensively during the last 20 years. Most important both types of models 
were found to be asymptotically free M . 

Starting from the early eighties so called RP 1 ^^ 1 models were discussed. 
The spins of these models are elements of the real projective space in N 
dimensions. This space can be thought of as a sphere S N_1 where opposite 
points are identified. Hence in perturbation theory the RP 1 ^^ 1 model is 
equivalent with the 0(A r )-invariant a- model. The fact however that the 
real projective space is not simply connected gives rise to topological defect 
structures similar to vortices in the 2D XY model. The questions discussed 
in the literature are whether these defects induce a phase transition at a 
finite coupling or whether these non-perturbative effects survive in the weak 
coupling limit. 

The lattice action of the RP N ^ model mostly discussed is 

S = ~P £ {S X Syf , (1) 

<xy> 

where < xy > is a pair of nearest neighbour points on the lattice and s is a 
unit vector in R N . An alternative way to identify — s x and s x is to introduce 
an Z2 gauge field 

S = —/3 Z <X y > S X Sy , (2) 

<xy> 

where z takes the values 1 or —1. 

Similar models have been introduced to describe orientational phase 
transitions in nematic liquid crystals ||]. These models have mainly been 
studied in three dimensions, where a weak first order phase transition is 
found (see ref. Q and references given in ||]). 

The numerical study of RP N ~ l models in 2D gave rise to much contro- 
versy. Some authors [||, ||, f?J find that their results are consistent with a 
phase transition at a finite coupling, while others doubt the existence of a 
phase transition but still see strong crossover effects between the strong and 
the weak coupling regime [Q, ||, ||] . 

Recently Caracciolo et al. (9|, [R], [ll|] argued that there is no phase 
transition in the RP 1 ^' 1 models. They claimed, based on a finite size scaling 
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analysis, that the RP 1 ^^ 1 models however have a weak coupling limit distinct 
from that of the 0(A r )-invariant <r-models. They claim even further, that a 
whole sequence of universality classes can be obtained from mixed models. 
In the following I will give evidence that rules out the scenario presented 
10, |ll]] . For a particular type of the action of the O(N) and the RP N ^ 
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model I will show that the models are exactly equivalent for sufficiently small 
coupling. I discuss the scaling properties of vortices of the RP N ^ model 
with standard action in the weak coupling regime. I give numerical results 
for the step scaling function introduced in ref. [12] that supports that the 
constraint model gives the same universal results as the standard action. I 
show that the differences in the finite size scaling curves for the RP 2 and 
the 0(3) model found in refs. [10, 11] can be explained as a boundary effect. 



2 The Constraint O(N) and RP 1 *- 1 Models 

Let me first define the models. The field variable s x is in both cases a 
unit vector in R N . In the case of the 0(-/V)-invariant model the Boltzmann 
weight of a configuration is equal to 1 if 

S X Sy > C (3) 

for all nearest neighbour pairs of sites < xy >. Else the Boltzmann weight 
is equal to 0. 

In the case of the RP N ~ l model — s x and s x are identified and the 
constraint on the field configuration is given by 

\s x s y \ > C (4) 

for all nearest neighbours < xy >. Equivalently one might introduce a gauge 
field 

z <xy> taking the values —1 or 1. 

> C . (5) 

In the following I shall show that the constraint 0(A r )-invariant model 
and the constraint RP N ~ 1 model are equivalent for C > cos(7r/4). Let us 
consider a lattice where all closed paths are contractible, i.e. all closed paths 
can be shrinked to an elementary plaquette by removing single plaquettes 
sequentially. A hyper-cubical square lattice with open boundary conditions 
is an example for such a lattice. 
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Consider the class of 2 , where V is the number of lattice points, con- 
figurations that arise from a given configuration s x by taking either +s x or 
— s x at each lattice point. Since 

\SxSy\ = \(-S X )Sy\ = \S X (-Sy)\ = \(-S X )(-Sy)\ (6) 

all configurations in such a class are either allowed or forbidden RP N ~ l 
configurations. Obviously a class of configurations that is forbidden under 
the RP N_1 constraint contains no configuration that is allowed under the 
0(N) constraint (with the same C). In the following I will demonstrate 
that for C > cos(7r/4) a class of configurations that is allowed under the 
RpN-i cons traint contains exactly 2 configurations allowed under the O(N) 
constraint, and therefore the partition functions are equal up to a trivial 
factor 2 y - 1 . 

Take one configuration out of an allowed class of RP N ~ l configurations. 
Pick one site x. Replace the spins on the other sites by 

8y' = 8y ]J sign(s u s;) . (7) 

<uv>& path(z,y) 

The result of this construction is independent of the paths chosen if 

JJ sign(s u s;) = 1 (8) 

<MD>Gclosed path 

for all closed paths on the lattice. For elementary loops consisting of four 
lattice-points this is the case for C > cos(-7r/4). All other paths can be 
successively built up out of elementary loops, since we have chosen a simply 
connected lattice topology. When adding an elementary loop the sign of a 
loop is conserved since the sign of the product of the new links in the path 
is the same as for the old links. Hence the sign of any closed path is 1 for 
C > cos(7r/4). 

The idea behind this proof has been discussed for an action similar to 
that in eq. (2) by Caselle and Gliozzi ||. However for that action the 
rigorous proof for a pure gauge in the weak coupling limit is missing. 

In Monte Carlo simulations one typically uses periodic boundary condi- 
tions, which leads to the lattice topology of a torus. Here loops exist that 
wind around the torus and hence cannot be contracted to an elementary 
loop. 
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In order to avoid configurations that are allowed under the RP con- 
straint but not for O(N) one has to require C > cos(tt/L) where L is the 
extension of the lattice in units of the lattice spacing. It is important to 
note that such boundary effects might well survive the continuum limit in 
a finite size scaling analysis. However this boundary effect can be repro- 
duced by proper boundary conditions imposed upon the 0(N) model. For 
C > cos(7r/4) a constraint RP N ~ 1 model on a periodic lattice is equiv- 
alent to a constraint 0(N) model with fluctuating boundary conditions. 
Fluctuating boundary conditions mean that in the partition function one 
sums over periodic as well as anti-periodic boundary conditions. In the 
case of anti-periodic boundary conditions one identifies s(0,y) = —s(L,y), 
s(L + 1, y) = —s(l, y) s(x, 0) = —s(x, L) and s(x, L + 1) = —s(x, 1). 



3 Scaling of the Vortex Density for the Standard 
Actions 

For the standard actions of the RP N ~ l model similar arguments apply. In 
the limit (3 — > oo the energy of a vortex should win against the entropy and 
vortices should play no role in the continuum limit of the theory. 

Let us identify a frustrated plaquette in eq. (2) with the center of a 
vortex. The classical solution of the s-field for a fixed gauge field with 
two frustrated plaquettes has an energy proportional to lnr where r is the 
distance in between these two frustrated plaquettes. Hence one can find a 
finite ro such that the energy is larger than 2&o + e, where &o is the leading 
coefficient in the perturbative /^-function. Therefore the density of vortex 
pairs with a distance larger than ro dies out faster than the square of the 
inverse correlation length. Hence they can not play a role in the continuum 
limit of the theory. 



4 Numerical Results for the Constraint Models 

In this section I show that the constraint O(N) model reproduces universal 
results of the 0(iV)-invariant a- model. Therefore I compute the step scaling 
function of ref. [^] for three different values of the running coupling and 
compare the result with that of ref. [12] obtained with the standard action. 
Furthermore I estimate the correlation length at C = cos(7r/4) using the 
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running coupling and also measure the correlation length for both the 0(N) 
and the RP N ~ 1 model for C < cos(7r/4) to check the importance of defects 
in the generation of the mass in the RP N ~ l case. 



The running coupling of ref. [12] is defined by 



9 2 = j-m(L)L , (9) 

where m(L) is the mass gap on a lattice with extension L in spatial direction. 
The /3-function for the running coupling g 2 is given by [|12]] 

„, N-2 a N-2 - (JV-l)(iV-2) „ , , 

The step scaling function a(s, u) is the discrete version of the ^-function. It 
gives the value of the coupling after change of L by a factor of s starting 
from a coupling u. 

The simulation was done using the evident modification of the single 
cluster algorithm A bond < xy > is called deleted if after the reflection 
of one of the spins s x or Sy the constraint s x Sy < C is still satisfied. 

A proof of ergodicity is given in the appendix. The simulation results 
listed in table 1 are based on about 10 7 single cluster updates. The correla- 
tion function was measured using the cluster-improved estimator [Gj| . The 
mass was extracted from the correlation function at distance L and 2L. 

Fitting the data of table 1 to an Ansatz 

S(2,n,a/L) = a(2,u) +c (a/L) 2 . (11) 

I obtain a(2, 1.0595) = 1.2589(10) from L/a > 16, a(2, 0.8166) = 0.9150(8) 
from L/a > 8 and a{2, 0.7383) = 0.8159(8) from L/a > 8. These re- 
sults can be compared with the step scaling function obtained in ref. |jl2: l 
er(2, 1.0595) = 1.2641(20), <r(2, 0.8166) = 0.9176(8) and er(2, 0.7383) = 
0.8166(9). The slight disagreement (about 2 standard deviations) might 
well be explained by deviations of the corrections to finite size scaling from 
the fit- Ansatz chosen. 

The exact prediction for the mass gap given by ]il 



mo , . 

= Z (12) 
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for N = 3 and the conversion factor for the A parameters 

-r'(i) 

47T 

given in ref. \12\ allows us to give an estimate for the infinite volume cor- 



A = ~^ A MS (13) 



relation length based on the measurement of the correlation length on a 
finite lattice. Taking the Monte Carlo result for the running coupling given 
in table 1 I obtain £ = 0.7 x 10 5 as estimate for the correlation length at 
C = 0.55 and £ = 0.6 x 10 9 as estimate at C = cos(7r/4), where the 0(3) 
and RP 2 constraint models become identical. 

In addition I performed some simulations for both the constraint 0(N) 
model and the constraint RP N_1 model at smaller C values such that the 
correlation length £ is much smaller than the lattice size L. Here I adopted 
the definitions used in refs. || [11], O . 

The correlation function in the vector-channel is defined by 

G v (x,y) = {s x s y ) . (14) 

Since naively this quantity vanishes identically under the symmetries of the 
rp n-i model 

one considers the tensor channel with the correlation function 

G t (x,y) = ((s x s y ) 2 )-^ . (15) 

Starting from these definitions of the correlation function one obtains the 
susceptibility 

and 

1 2ir 

F = yYl ™s( — k(x - y))) G(x, y) (17) 
x,y 

with k = (1,0) or k = (0, 1). 

The second moment correlation length is now defined as 

^ " 2sin(vr/L) ' (18) 

In table 2 some results for the constraint 0(3) model are given. It is 
remarkable that already for C = the correlation length is larger than ten. 
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The ratio of the correlation length in the vector and the tensor channel is 
about £u/£t = 3.3(1). 

In table 3 my results for the constraint RP 2 model are summarized. At 
C = 0.55 there is a factor of about 10 3 in between the correlation lengths of 
the 0(3) and the RP 2 model. This means that it is practically impossible 
to see the true asymptotic behaviour of the constraint RP 2 model in a 
computer simulation. 



5 Finite Size Scaling and Universality 

In this section I shall demonstrate that the difference in the finite size scal- 
ing curves observed in ||] and |1C] can be explained in part by the boundary 
effect discussed above. I simulated the 0(3)-invariant model with the stan- 
dard action on a square lattice using fluctuating boundary conditions in 
both lattice-directions. For the updates of the boundary conditions I used 
the boundary flip algorithm proposed in ref. |l6| for the Ising model and 
generalized to O(N) vector models in ref. [jT^| . 

I performed runs at (3 = 1.4, (3 = 1.5 and (3 = 1.6. The true correlation 
lengths for these (3 values are f = 6.90(1), £ = 11.09(2) and £ = 19.07(6), re- 
spectively flq| . I used lattice sizes ranging from L = 6 to L = 128. Through- 
out I performed 100000 measurements. I performed a measurement after one 
boundary-flip update for each direction and roughly cluster-size/ lattice- 
size standard single cluster updates. In figure 1 the dimensionless quantity 
Xt(2L) I Xt{L) is plotted as a function of £/L. I give the results for fluc- 
tuating boundary conditions (circles) and for comparison the results with 
periodic boundary conditions (diamonds). In order to compare with fig. 2 
of ref. |ll| one has to take the factor £ v /£t = 3.3(1) into account. My result 
for periodic boundary conditions is consistent with that given in fig. 2b of 



ref. [11]. The fluctuating boundary conditions remove the characteristic dip 
visible in the finite size scaling curve for periodic boundary conditions. The 
finite size scaling curve for periodic boundary conditions looks qualitatively 



much like that of fig. 2a of ref. 11] (RP 2 like models). However the slope 
of the curve in fig. 2a of ref. [O] is much steeper for large £/L than that of 
fig. 1. This shows that the results given in fig. 2a of ref. [|llj are effected 
by strong corrections to scaling due to vortices. 
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6 Is There a Phase Transition? 



In order to understand the phase-structure of the RP models one might, 
in analogy with the KT scenario of the XY model |18| ], discuss the RG- 
flow of the models in a 2 dimensional parameter space. In addition to the 
coupling g 2 one might introduce a coupling parameter \i for the plaquette- 
term, controlling the density of vortices. 

S= ^ Z xy S x Sy + 111 fl^Zp , (19) 

9 <xy> p 

where z p = U<xy>e P z xy 

I will make no attempt here to derive the RG flow-equations. However 
certain qualitative features and their consequences seem to be evident: 

a) For (g 2 ,0) the standard /3-function of the O(N) model is recovered. 

b) Vortices cause disorder. Therefore a non-vanishing fj, should amount 
to a positive contribution in the derivative of g 2 with respect to the logarithm 
of the cutoff scale and hence accelerate the flow towards strong coupling. 

Statement a) rules out that a possible phase transition in RP N ~ 1 is KT 
like, since the fixed-point of the if T-transition is purely Gaussian. Further- 
more statement b) rules out any fixed-point that might occur at a finite 
p. 

Still we have to explain why Monte Carlo simulations and strong coupling 
expansions seem to be in favour of a phase transition. It seems plausible 
that in analogy with the if T-flow equations fi is irrelevant for small coupling 
g 2 but becomes relevant above some threshold value g 2 . That means above 
g 2 the RG-trajectories are driven off from the renormalized trajectory of the 
0(A r )-invariant model. 

7 Conclusions 

I have proven that the constraint O(N) and constraint RP 1 ^' 1 model be- 
come equivalent for C > cos(7r/4). Using the renormalized coupling g 2 = 
m(L)L I estimated the correlation length at C = cos(-7r/4) to be about 
£ = 0.6 • 10 9 for both models with N = 3. For C-values being smaller, such 
that £ << 1000, the models display huge differences. This means that the 
asymptotic behaviour of the constraint RP 2 practically can not be observed 
in a computer simulation. I argue that a similar scenario holds for models 
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with a standard action. As (3 — > oo vortices in the RP N ~ 1 model vanish and 
the RP N ~ 1 becomes equivalent to an O(N) model by the virtue of a gauge 
fixing. 

On lattices with periodic boundary conditions one has to notice that 
paths winding around the lattice are not contractible. The effect of such 
loops in the RP N_1 model amount to fluctuating boundary conditions in the 
equivalent O(N) model. I demonstrated numerically that this fact partially 
explains the differences found in the finite size scaling curves for the 0(3) 
and RP 2 models observed in refs. [10, O]. 



8 Note added 



My conclusions are confirmed by the work of F. Niedermayer, P. Weisz and 
D.-S. Shin [19], which I found today on the hep-lat bulletin-board. 



9 Appendix 

In the following I prove that the single cluster algorithm applied to the 
constraint 0(A r )-invariant model is ergodic. 

It is sufficient to show that any allowed configuration can be transformed 
in a finite number of cluster-update steps to the configuration s = (1, 0, 0) 
for all sites. 

Let us consider a N = 2 (XY) model with a bond dependent constraint 
C <xy> . Assume that the spins are distributed in an angle range [0, a&] with 
respect to the 1-axis. (The largest range to start with is [0, 2n].) 

Take a reflection axis which has an angle a&/2 with the 1-axis. Per 
construction none of the sites x with (j) x > a^/2 is connected via a frozen 
bond with a site y with (fi y < at/2. Hence all spins can be moved into the 
range [0, + 1] with otk+i = a-k/2 using a finite number (smaller or equal 
the number of sites) of cluster updates. Iterating this process, in a finite 
number of steps all spins can be put into the range [0, minarccos(C <a; y>)]. 
Now take for each site a reflection axis with a x = 4> x /2. Per construction 
all these clusters are single site clusters. 

We hence constructed a sequence of a finite number of cluster updates 
that transforms an arbitrary configuration to the s = (1,0) for all sites 
configuration. 
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For general N this procedure can be iterated. Consider the N and (iV— 
l) th component as an embedded XY model. Remove the N th component. 
Go ahead untill s = (1,0, ...,0) for all sites is reached. 
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Figure captions 



Figure 1: The dimensionless quantity Xt(2L) /xt{L) is given as a function 
of £ V /L for the 0(3) model. The data points with the circles are obtained 
with fluctuation boundary conditions while those with diamonds are ob- 
tained with periodic boundary conditions. 
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Table 1: The renormalized coupling g 2 from the constraint 0(3) model. C 
gives the value of the constraint. L/a and V ja are the lattice extensions in 
spacial direction. 
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Table 2: The second moment correlation length in the vector (£„) and tensor 
(£t) channel for the constraint 0(3)-invariant vector model for various values 
of the constraint C. 
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Table 3: The second moment correlation length in the tensor channel £j for 
the constraint RP 2 model for various values of the constraint C . 
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